ELEMENTARY INVERSION OF RIESZ POTENTIALS 
AND RADON-JOHN TRANSFORMS 



B. RUBIN 

Abstract. New simple proofs are given to some elementary ap- 
proximate and explicit inversion formulae for Riesz potentials. The 
results are applied to reconstruction of functions from their inte- 
grals over Euclidean planes in integral geometry. 



1. Introduction 

The Riesz potential of order a of a sufficiently good function / on 
M" is defined by 

Rea > 0, a n,n + 2,n + 4, . . . . 

This operator can be regarded (in a certain sense) as a negative power 
of "minus-Laplacian" , namely, 

2 d'^ (9^ 9^ 

^ ^ ' ^ dxf ~^ dxl ~^ ' ' ' ~^ dxf^' ^^ '^^ 

If f E L^(M'^), then the integral (11.11) converges a.e. provided that 
1 < p < n/Rea. This condition is sharp. Explicit and approxi- 
mate inversion formulae for /° are of great importance. We refer to 
[3 [121 ED 1231 ESI EZl ED] for the basic properties and applications of 
Riesz potentials. In this article we restrict to real a G (0,n). This case 
reflects basic features and is sufficient for integral geometric applica- 
tions in Section HI 

Our main working tools are a suitably chosen auxiliary function w{x) 
(we call it a reconstructing function) and its scaled version Wt{x) = 
t~"'w{x/t), t > 0. In the modern literature such a function is called a 
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wavelet and a convolution (/ * Wt){x) is called the continuous wavelet 
transform of /; cf. |21 E]- We choose 

w{x) = (-A)"[(l + w{x) = {-Aw){x), (1.3) 

m G N, and keep this notation throughout the paper. 

The main results of the article are presented by Theorem A below 
and Theorem B in Section HI The latter gives an example of applica- 
tion of Theorem A to elementary inversion of the fc-plane Radon- John 
transform in integral geometry and falls into the general scope of the 
convolution-backprojection method (cf. 122] )• Radon transforms 
and their generalizations are the basic tools in tomography and numer- 
ous related areas of pure and applied mathematics; see [H El El El El El 
Uni m in 11511201 [11] and references therein. 

1 < p < n/a. Then 

Ca,m = InC^m - a), (1.4) 
da,m = {2m-a) Ca,m, (1-5) 

where /q°°(. . .) = lim ■ ■) and 7n(-) has the same meaning as in 

( li. jj) . The limit in both formulae exists in the L^-norm and in the a.e. 
sense. If, moreover, f G Co(M"), then it exists in the sup-norm. 

We recall that Co(M") denotes the space of continuous functions on 
M" vanishing at infinity. 

Some comments are in order. Statements of this kind are not new. 
Regarding (11. 4p , we observe that approximate inversion of operators of 
the potential type was initiated by Zavolzhenskii, Nogin and Samko [321 
[T71 1211 [251 126] . Conceptually this approach is close to the convolution- 
backprojection method in tomography |T5]. An elegant formula (II. 4p 
is due to Samko [26l p. 325], [241 EH]; see also [19]. Below we give a 
new simple proof of it. 

Regarding (II. 5p . wavelet-like inversion formulae for Riesz potentials 
are also well-known; see, e.g., [211 Section 17]. Our aim is to show that 
elementary reconstructing functions (II. 3p work perfectly in the wavelet 
inversion scheme and the relevant justification is much simpler than 
in the general theory; cf. [21]. The case a = in (II. 5p . when J" is 
substituted by the identity operator, represents a variant of Calderon's 
reproducing formula [21 E]. 

The validity of formulae (II. 4p and (II. 5p can be formally explained in 
the language of the Fourier transform. We have 

lim (r/ * t-'^wtno = Ho mr'^wito = /(o (i.e) 



Theorem A. Let 2m> a, f e U'{^ 
c^,mf = lim(r/*t-"w;i) 
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provided that the hmit Cw = hm exists. Moreover, since w 

in (11. 3p is a radial function, we easily get 



^ I''f*wt^^\^, , , 1 [ w{y)dy 



(0=^^/(0, d^ = ^ y^Sdl' (1-7) 

On-l JM" |y| 



where o"„_i =27r"/^/r(n/2) is the surface area of the unit sphere |x| = 1. 

This simple argument does not work for arbitrary L^-functions. To 
cover this case, we develop another technique and give the final answer 
without using the Fourier transform. 



2. Preliminaries 

We recall some known facts. The Fourier transform of a function 
f eL^ = Li(R") is defined by 

/(0=/ /(x)e"-«dx, x-e = Xiei + ...+:r„en. (2.1) 

Let S = ^(R") be the Schwartz space of all C°° functions which, to- 
gether with derivatives of all orders, vanish at infinity faster than any 
inverse power of |x| = {xl + . . . + x'^Y^'^. We endow S'(R"') with a 
standard topology generated by the sequence of norms 

||v?|U = max(l + |x|r y m = 0,1,2,... . 

The following spaces adapted to Riesz potentials were introduced by 
Semyanistyi [28]; further generalizations due to Lizorkin and Samko 
can be found in |l3l |26] . 

Let = \E'(R'^) be the subspace of S, consisting of functions ipiC.) 
vanishing at ^ = with all derivatives, and let $ = $(R") be the 
Fourier image of We equip $ with the induced topology of S. Then 
$ becomes a topological vector space which is isomorphic to \E' under 
the action of the Fourier transform. We denote by $' the space of 
distributions over $. 

The main reason for introducing the spaces and $ is that \E' is 
invariant under multiplication by |^|" for any a G C and therefore, the 
Riesz potential J" is an automorphism of $ and 

Proposition 2.1. [21, p. 21 ] Let f E L'' , I < r < oo, and g E L'p , 

1 < p < oo. If f = g in the -sense, then f = g almost everywhere. 
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3. Proof of Theorem A 

The following auxiliary propositions form the heart of the proof. 

Lemma 3.1. Let weL^(R"') be such that h = I°'w has an integrable 
decreasing radial majorant. If f eLp{M."'), l<p<n/a, then 

\im {r f *t~°'wt) = Cnj{a) f, c^(a) = / h{x)dx, (3.1) 

where the limit exists in the U'-norm and in the a.e. sense. If, more- 
over, f G Co(M"), then the limit exists in the sup-norm. 

Proof. We have = t-^hix/t) = ht{x). Then 

rf * t-'^wt = f * t-'^rwt = f*ht, 

which yields an approximate identity [SOj p. 62]. The above application 
of Fubini's theorem is permissible because (/"|/|) * |w| < oo a.e. □ 

There are many ways to choose a reconstructing function w. Usually 
w or /"w or both are expressed analytically in a pretty complicated 
way. As we shall see below, an advantage of the choice (11. 3p is that both 
w and h = I'^w are elementary and a constant Cw{(y) can be different 
from zero. 

Lemma 3.2. Let w{x) = (-A)™[(1 + |a;|2)m-(n+a)/2] ^ If a > 0, then 
w&L^iW^). If, moreover, a < min(n, 2m) and h = I"w, then 

_^r{{n + 2m-a)/2) 



r((n+a — 2m)/2) 

so that h e Li(R") and 

22'"-°7r"/2r(m-a/2) 

Ca,m= h{x)dx= — — — = 7„(2m-a); 

Jun T{{n + a-2m)/2) 

cf. the normalizing constant in ( li.ij) . // a is not an integer, then 
Ca,m 7^ for any m > a/2. // a G {1, 2, . . . , n — 1}, then Ca,m 7^ 
provided, e.g., that m = [(a + 2)/2]. 

Proof. The first statement can be easily checked by differentiation. For 
instance, we can write A in polar coordinates to get 

w(x) = (L"'^Ij)(\x\^), L = 4ri-"/2A^«/2^^ V^(r) = (l+r)'"-("+")/2_ 

dr dr 

This gives w{x) = 0((1 + IxH-^^+^^Z^) ^ L^(M"). To prove the second 
statement we invoke the Bessel-McDonald kernel 

^"^^^ " |^|(n-")/2 ' = W2r(a/2)' ^^'^^ 
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where Ku{-) denotes the modified Bessel function (or the McDonald 
function) of order u with the property 

K,{z) = K^,{z)- (3.3) 

see, e.g., [3lj. It is known that 

'|^|(a-n-l)/2g-|€| if|^|>l, 

ler-" if lel < 1, a < n, 

1 if 1^1 < 1, a > n, 

[l + log(l/|^|) if |e|<l, a = n. 



(3.4) 



where Cq, is cL continuous function of a; see, e.g., [211 P- 257], [HI p. 
285]. Furthermore, for any a > the Fourier transform of (1 + |a; 
in the sense of distributions is computed by the formula 



2\-a/2 



((1 + |a;^-"/^0(a;)) = (27r)-"(G,(0, 0(0), e ^(M"), (3.5) 
[m Section 8.1]. Let us prove that 

{rw, <P) = a„,„ ((1 + |a;|2)(-2— n)/2^ ^(3,)) (3.5) 

for any test function (p in the space $(M"). Once this is done, the 
required pointwise equality will follow owing to Proposition 12.11 We 
have 

{rw, 0) = {W, r(P) = ((1 + |a;|2){2m-n~a)/2^ {-Ay''{I"(f)){x)) 

= (27r)-"(G„+„_2™(O,iel'"'"0(O) 
= (2vr)-"A„+„_,^ 0(0 



Using 03. 3p and fl3.5p . we continue: 



(rwA) = ^^4-^^^^^^(G2™+n-„(O,0(O) 

= a„,^((l+|x^("-^"-")/^0(x)). 

Evaluation of c^,™ is straightforward. To complete the proof, we note 
that Ca^m 7^ if (n + a — 2m)/2 7^ 0, —1, —2, .... The latter is guaran- 
teed under the afore-mentioned choice of m. □ 

Proof of Theorem A. The first statement follows immediately from 
Lemmas 13.11 and 13.21 To prove the second statement, we first note 
that w = —Aw belongs to L^. Hence, the convolution J"/ * t'^'wt is 
well-defined in the Lebesgue sense and can be written as f * ht with 
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h = I°'w = /"(— A)w. Furthermore, for any test function G $(I 
dSSD yields 



Ch,<P) = (r(-AV,0) = -(/-«;, A 0) 

= -a„,™(A[(l + |x|^)("-")/2-'"],0). 

Since A [(1 + |a;|2)(a-n)/2-m] = o{{l + |a;|2)("-")/2-"^-i) G then, by 
Proposition I2.H we have a pointwise equahty 

Denote 



Then 



Je = y ——dt = f*%l)e, iJe{x) = J —^dt- 

Setting h{x) = /io(|a;p), we get 



ipsix) 



2 X " /n ^ 



.|2 ~ 



This is a scaled version of the function = 2 ^\x\ " /q ho{r)r"^^'^ ^dr. 
Observing that 

ho{r) = -a^,m L [(1 + r)("-")/2-], L = 4ri-"/2 



we have 



= -aa,m (a - n - 2m) (1 + |a;|2)("-n)/2-m-i_ 
This gives lim = damfix), where 

22™-"+i7r"/2r(m + l-a/2) , , , 

da,m = ^{x)dx= — — = (2m-a) 7„(2m-a), 

Jr„ r((n + a-2m)/2) 

and the limit is understood in the required sense. 
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4. Inversion of the Radon-John Transform 

We recall basic definitions. More information can be found in [5l [71 
Let Qn,k and Gn,k be the affine Grassmann manifold of all non- 
oriented /c-dimensional planes r in R" and the ordinary Grassmann 
manifold of fc-dimensional linear subspaces ( of M", respectively. Each 
fc-plane r G Qn,k is parameterized as r = (C, u), where ( G Gn,k and 
M G C"*" (the orthogonal complement of ( in M"). The k-plane Radon- 
John transform of a function / on M" is defined by 

(i?fc/)(r) ^ (i?fc/)(C, u) = J^f{y + u)dy, (4.1) 

where dy is the induced Lebesque measure on the subspace C G Gn,k- 
This transform assigns to a function / a collection of integrals of / over 
all fc-planes in M". The corresponding dual transform of a function ip 
on Qn,k is defined as the mean value of ip (r) over all A;-planes r through 
X G M": 

{Rl^)ix)= [ ^iaCo + x)da, x G M". (4.2) 

Here (q G G„,fc is an arbitrary fixed /c— plane through the origin. If 
/ G Lp(M"), then R^f is finite a.e. on Qn,k if and only ii 1 < p < n/k 



A variety of inversion procedures are known for Rkf', see 
One of the most important algorithms relies on the Fuglede formula [1] , 

Rk^kf = dk,nl''f, dk,n = {'2n)^ an-k-l/(^n-l, (4.3) 

with the Riesz potential I^f on the right-hand side. Hence, Theorem 
A yields the following result. We denote 

M;(x) = (-A)'"[(l+|xn™-("+'=)/2], m=[{k+2)/2], w{x) = {-Aw){xy, 

4m^(n+fc)/2r(^/2) r(m - fc/2) 
- Tiin-k)/2)Tiin+k)/2~my " ^^"^'^^ 

Theorem B. If = Rkf, f G Lp{W), l<p< n/k, then 

Xkf = \iMRlip*t-''wt), 6k f= -^^^dt, (4.4) 

where J^{. ■ •) = lira JJ^(. . .). The limit in both formulae exists in the 

L^-norm and in the a.e. sense. If, moreover, f G Co(M"), then it exists 
in the sup-norm. 
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